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Variational formulation of the motion of an ideal fluid on the basis of gauge principle 
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On the basis of gauge principle in the field theory, a new variational formulation is presented for 
flows of an ideal fluid. The fiuid is defined thermodynamically by mass density and entropy density, 
and its fiow fields are characterized by symmetries of translation and rotation. A structure of 
rotation symmetry is equipped with a Lagrangian Aa including vorticity, in addition to Lagrangians 
of translation symmetry. From the action principle, Euler's equation of motion is derived. In 
addition, the equations of continuity and entropy are derived from the variations. Equations of 
conserved currents are deduced as the Noether theorem in the space of Lagrangian coordinate a. It 
is shown that, with the translation symmetry alone, there is freedom in the transformation between 
the Lagrangian a-space and Eulerian a;-space. The Lagrangian Aa provides non-trivial topology 
of vorticity field and yields a source term of the helicity. The vorticity equation is derived as an 
equation of the gauge field. Present formulation provides a basis on which the transformation 
between the a space and the x space is determined uniquely. 
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I. INTRODUCTION 

In the historical paper ' General laws of the motion of 
fluids^ [1], Leonhard Euler verified that his equation of 
motion can describe rotational flows. The same theme is 
investigated in this paper under modern view. Fluid me- 
chanics is understood to be a field theory in Newtonian 
mechanics that has Galilean symmetry. It is covariant 
under transformations of the Galilei group. The gauge 
principle ([19], [20], [21]) requires a physical system under 
investigation to have a symmetry, i. e. a gauge invariance 
with respect to a certain group of transformations. Fol- 
lowing this principle, the gauge symmetry of flow fields 
is studied in [2] and [3] with respect to both translational 
and rotational transformations. The formulation started 
from a Galilei- invariant Lagrangian of a system of point 
masses which is known to have global gauge symmetries 
with respect to both translation and rotation [4]. It was 
then extended to flows of a fluid, a continuous material 
characterized with mass density and entropy density. In 
addition to the global symmetry, local gauge invariance 
of a Lagrangian is required for such a continuous field. 
Symmetries imply conservation laws. Equations of con- 
served currents are deduced as the Noether theorem. 

Thus, the convective derivative of fluid mechanics, i.e. 
the Lagrange derivative, is identified as the covariant 
derivative, which is a building block in the framework 
of gauge theory. Based on this, appropriate Lagrangians 
are defined for motion of an ideal fluid. Euler's equation 
of motion is derived from the action principle. In most 
traditional formulations, the continuity equation and en- 



tropy equation are given as constraints for the variations, 
while in this new formulation those equations were de- 
rived from the action principle. In the previous study 
([5], [6]) of rotational symmetry of the velocity field v{x), 
it is found that the vorticity u) — V xv is the gauge field 
associated with the rotational symmetry of velocity. 

A new structure of the rotational symmetry was given 
in [3] by the following Lagrangian: 
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where A is a vector potential and CwA = dtA-\-v''dkA-\- 
AkVv'^. This is derived from a representation character- 
istic of a topological term known in the gauge theory. This 
yields non-vanishing rotational component of the veloc- 
ity field, and provides a source term of helicity. This 
is closely related to the Chern-Simons term, describing 
non-trivial topology of vorticity field, i.e. mutual linking 
of vorticity lines. The vorticity equation is derived as an 
equation for the gauge field. 

With regard to the variational formulation of fiuid 
flows, the papers [7] and [8] are among the earliest to have 
influenced current formulations. Their variations are car- 
ried out in two ways: i.e. a Lagrangian approach and 
an Eulerian approach. In both approaches, the equation 
of continuity and the condition of isentropy are added 
as constraint conditions on the variations by means of 
Lagrange multipliers. The Lagrangian approach is also 
taken by [9]. In this relativistic formulation those equa- 
tion are derived from the equations of current conser- 
vation. Several action principles to describe relativistic 
fluid dynamics have appeared in the past (see [9, §4.2] 
for some list of them). 

In the Lagrangian approach, the Euler-Lagrange equa- 
tion results in an equation equivalent to Euler's equation 
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of motion in which the acceleration term is represented as 
the second time derivative of position coordinates of the 
Lagrangian representation. In this formulation, however, 
there is a certain degree of freedom in the relation be- 
tween the Lagrangian particle coordinates and Eulerian 
space coordinates. Namely, the relation between them 
is determined only up to an unknown rotation. In the 
second approach referred to as the Eulerian description, 
the action principle of an ideal fluid results in potential 
flows with vanishing helicity, if the fluid is homentropic 
[2]. However, as noted in the beginning, it should be pos- 
sible to have rotational flows even in such a homentropic 
fluid. Gauge theory for fluid flows provides a crucial key 
to resolve these issues. It was also shown in [2] that a 
general solution in the translational symmetry alone is 
equivalent to the classical Clebsch solution [10]. A new 
formulation on the basis of the Clebsch parametrization 
is carried out in [11] and [12] aiming at its extension to 
supersymmetric and non-Abelian fluid mechanics. 

It is interesting to note the gauge invariances known 
in the theory of electromagnetism and fluid flows. There 
is an invariance under a gauge transformation of electro- 
magnetic potentials consisting of a scalar potential cf) and 
a vector potential A. An analogous invariance is pointed 
out in [2] for a gauge transformation of a velocity po- 
tential (j> of irrotational flows of an ideal fluid, where the 
velocity is represented as = Wcj). It is shown in [7] 
(cited in [3]) that gauge invariance is not restricted to 
the potential flows, but also there is known an invariance 
in the rotational flow of Clebsch representation. 



II. EQUATIONS IN o-SPACE 
A. Lagrangian 

Let us consider a variational formulation with a La- 
grangian represented with the particle coordinate a = 
{a^,a'^,a^) = (a, 6, c) {i.e. Lagrangian coordinates) . In- 
dependent variables are denoted with a'^ where ii or greek 
letter suffix take = 0,1,2,3 with a° the time variable 



written also as t {= t): 
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a'^). Correspond- 



ing physical space coordinate x = (x. y, z) {Eulerian co- 
ordinates) are written also as = {t 3). The 
letter r is used (instead of t) in combination with the par- 
ticle coordinates a''. Physical space position of a particle 
a is expressed by X'=(a'^') = X'=(r, a), or X'^ = {X, Y, Z). 
Its velocity is given by v'^ = drX'^, also written as X!^. 

The Lagrangian coordinates (a, b, c) are defined such 
that an infinitsimal three-element d^a = dad6dc denotes 
a mass element dm of an infinitesimal volume d^x = 
dxdydz of the cc-space. The mass elelment dm should 
be invariant during the motion: 

dr{dm) = dr{d^a)=0. (1) 

The mass-density p is defined by the equation d^a = 



transformation X'^ = X*^(a') from a-space to X-space 
{k, I = 1,2, 3), we have 



d{X\X^,X^) _ d{X, Y, Z) 
5(a^,a2,a3) d{a,b,c) 



(2) 



In an ideal fluid, there is no dissipation of kinetic energy 
into heat, by definition. According to thermodynamics 
for the entropy s (per unit mass) and temperature T , we 
have T(5s = if there is no heat production. Namely the 
entropy s does not depend on r. Then, the change of 
internal energy e (per unit mass) is related to the density 
change 5p alone by 
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6h- 



1 



6p, (3) 



where p is the fluid pressure, and h = e+p/ p the enthalpy, 
and ( • )s denotes s being fixed. However, the entropy s 
may not be uniform and may depend on a by initial 
condition. Hence, s = s{a), or equivqalently, 



drS = 0. 

Total Lagrangian is defined by 



At = / i X'; X'; d'^a - / e{p, s) d^a , 

JMa JMa 



(4) 



(5) 



[3] , where Ma is a space of fluid under investigation, and 
X!^ = Xq = v'^ is the velocity. The internal energy e(p, s) 
of the second term depends on p (which in turn depends 
on X^ = dX^jda^ by (2)) and the entropy s{a). 
An action 7 is defined by the integral: I = J Axdr: 

I = j L{X^)d\ d'^a = dTd^a, (6) 
i(X^) = iXo^Xo^-e(Xf,a'=). (7) 



B. Noether's theorem 

Euler-Lagrange equation associated with the La- 
grangian (7) is given by 



d r dL 



dL 
'dX^ 



dL 
AdX^ 



dj 



dL 
'dX^ 



= 0. 



Energy-momentum tensor is defined by 



(8) 



(9) 



[7], where k = 1,2,3. As long as (8) is satisfied to- 
gether with an assumption of r-independence of L {i.e. 

d-T-L = 0), it can be verified [3] that we have a conserva- 
tion equation d^T^ = (where 3^ = d/da'^). This is the 
Noether theorem ([13], [19]). 

For /i 7^ (.T^ = a), the conservation law Qi/T^ = 
reduces to the momentum equaitons: 



pd X. With using a Jacobian determinant J of the drVa + daF = {Va = XaXr + Yo,Yr + ZaZr), (10) 
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[7], where F = —^v +h. Two other equations are ob- 
tained with a replaced by cychc permutaion of {a,b,c). 
Integrating this with respect to r between and t, we 
find the Weber's transformation [14, Art. 15]: 

F„(t) = X^X, + Y^Y, + Z^Z, = V^{0) - d^x, (11) 



/ FAt= [ (-it;2 + /i)dT. 
Jo Jo 



The Va of (10) is a transformed velocity in the a-space 
(Sec.V A). Its time evolution is given by (11) for a given 
initial values of V^(0, a) and h{0, a) at a = x. 
With /z = 0, we have the energy equation: 



drH + da 
+ dc 



d{X,Y,Z) 



d{T, b, c) 
d{X,Y,Z) 



d{a,b,T) 



+ db 
= 0. 



d{X,Y,Z) 



d{a, T, c) 



(12) 



where H = ^v"^ + e. The equation (10) reduces to the 



equation for the acceleration Aa{T, a): 



XryX-r 



■ daP, (13) 



which is known as the Lagrangian form of equation of 
motion [14, Art. 13]. This can be transformed to 



1 
p 



dxP- 



da dp 
dx da 



(14) 



[3]. Since Xtt is the x-acccclcration of the particle a. 
this is the form equivalent to the x-component of Euler's 
equation of motion (25). The y and z components can 
be obtained analogously. 



C. Arbitrciriness in the transformation 

There is an arbitrariness in the transformation from 

the a-spacc to the a;-space with respect to the equa- 
tion (13). Its middle-side expression is a form of scalar 
product of two vectors in the a;-space: the particle 
accelcrataion (X^^, 1^^, Z^^) and the direction vector 
[Xa, Ya,Za) of the a-axis in the a-space. 

Putting it in a different way, the equation (13) is in- 
variant with respect to orthogonal rotational transforma- 
tions of a displacement vector AX = {AX, AY, AZ) of 
a particle in the a;-space. In fact, suppose that a vector 
AX satisfies the equation (13). Then, another vector 
AX = R AX obtained by an orthogonal transformation 
R satisfies the same equation, since any orthogonal ma- 
trix satisfies RR^ = I (unit mtrix) where R^ denotes 
the transposed matrix of R. So that the vector AX is 
not uniquely determined. The same freedom can be said 
to the velocity Va{T, a) of (11) as well. 

These imply that a certain machinery must be 
equipped in order to fix this arbitrariness within the 
framework of rotational symmetry. This will be consid- 
ered later. Note that the density p is not changed by the 
orthogonal transformation. 



III. EQUATIONS IN a;-SPACE 
A. Action in Eulerian representation 

Eulerian description is represented by the indepen- 
dent variables {t, x, y, z). Local gauge symmetries of fluid 
flows are investigated in detail in [2] , [3] . The time deriva- 
tive dr is equivalent to the convective derivative D(: 

9x = Dt, T)t = dt + udx + vdy + wdz = dt+v-V . (15) 

The operator T>t is verified to be gauge-invariant. The 
velocity field v{x, t) is defined by the particle velocity: 

v{x,t)=dTX = T)tX. (16) 

The acceleration field A{x, t) is also defined by 

A{x, t) = dT^X = T)tv = {dt + v''dk)v. (17) 

As noted previously, the mass d^a(a) and the entropy 

s = s{a) satisfy (1) and (4). In view of these properties, 
we can define the following two Lagrangians: 



Ld, = — I dT(l>^a, = — I adri^^a 

Jm Jm 



(18) 



where (/>(a, r) and il^{a, r) are scalar fields associated with 
mass and entropy, respectively. By adding and to 
At of (5), the total Lagrangian is given by 



j dT(l>d^a- j sdT^ d^a. 



(19) 



The action is defined by / = J^^ Arp*dr, where the inte- 
gral ^0 = / dr / dT(j)d^a can be integrated with respect 
to r and expressed as J[(j)]d^a, where [cj)] = 4>\t2 — ^In is 
the difference of (t> at the end times T2 and ti and hence 
independent of r G (ti,t2). Likewise, the last integral 
can be expressed as ~ J['ip]sd^a, because s is inde- 
pendent of T. This means that the gauge potentials (f> 
and ijj do not appear in the equation of motion obtained 
through variations of the action 7 for r S (ti,T2). 

However, it becomes soon clear that these are non- 
trivial in the expressions of the a;-space, because they 
are rewritten as = — Jj^ p Dt(^ d^a;, and = 
— ps 'Dt'tp d^x by using the relations d^a = pd^x and 
dT = T)f 

In the a;-space, the total Lagrangian can be written as 
Arp* = Jj^C{v,p,s,(j),'il)) d^x, where 

/:=^pv''v''- pe{p, s) -pDt<l>- psDt V (20) 

[22] . This is proposed as a possible form of Lagrangian in 

the a;-spacc (but an additional term will be added later). 
The action is defined by / = / C{v, p, s, 4>, ^) d'^x, where 
d^x = dtd^x. However, the action principle results in the 
potential fiow represented by i; = gTa,d{(j)+so'il>) when the 
fluid has a uniform entropy sq (see [2]). 
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B. Outcomes of variations 

We require invariance of the action / with respect to 

variantions. First, consider the following infinitesimal 
transformation: x'{x,t) = x + ^{x,t). The volume ele- 
ment d^x is changed to d^a;' = (1 + dk^'')d^x, up to the 
first order terms. Hence the variation of volume is given 
by A(d'^a;) = dkS.'^d^x, while the variations of density, 
velocity and entropy are Ap = —p dk£,'^, Av = Dt^, and 
As = 0. Under these together with (1) and (4) (with 
keeping </> and tp fixed), the variation of / is given by 



A7 



-Av + -Ap+-As + Ldke 



This is required to vanish for arbitrary variation of 
which results in the Euler-Lagrange equation: 

Similarly, invariance of / with respect to arbitrary vari- 
ations of (p and tp (denoted by Acp and At/j) leads to 

A(j) : dtp + V-{pv) = (continuity equation), (22) 
AV- : dtips) + V • (psv) = . (23) 



C. Noether's theorem in Eulerian representation 

Associated with (21), one can define the momentum 
density rrik and momentum-flux tensor by 



ruk 



Ml 



Prom (7), we obtain mk = pvk and Ml = pVkv'' +pSl., 
where Vk = v^' in the present Eucledian space. The equa- 
tion (21) can be written in the form of momentum con- 
servation, dt{pv'') + di[pv''v'^) + dkP = (dk = d/dx^). 
Using (22) , this equation can be reduced to the following 
Euler's equation of motion: 



--dkP {=-dkh). (25) 



The equation (14) is equivalent to this equation. 

The energy equation (12) can be transformed to the 
following equation of energy conservation: 

dt[p{^v^ + e)] +dk[pv''{\v^ + h)] =0. 



IV. ROTATION SYMMETRY 

A topological structure of vorticity field is now consid- 
ered with respect to the rotational symmetry. Related 
gauge group is the rotation group 50(3). An infinitesi- 
mal rotation is described by the Lie algebra so(3) of three 
dimensions, which is non-Abelian. 



From the study of the rotational gauge transformation 
[3], it is found that the covariant derivative Vf, velocity 
V and accerelation A are represented as 

Vt = dt + {vV), (26) 

v = VtX = {dt + {v-V))x, (27) 

A = \7tv = dtV + {v-V)v (28) 

VtV = dtv + gra,d{^ v^) + u: X V . (29) 

It is verified that the last expression of VtV = dtv + 
V{^v'^) + u; X V not only satisfies the rotational gauge- 
invariance, but also expresses that a; is the gauge field 
of the rotational symmetry. In addition, it satisfies the 
covariancc requirement with respect to Galilean trans- 
formation from one reference frame {t, x, v) to another 
(t*, X*, V*) moving with a uniform relative velocity U, 
where t^, = t, x^ = x — Ut and f , = v — U. Namely, we 
have the covariance VtV = (Vtt))*. 



V. LAGRANGIAN ASSOCIATED WITH 
ROTATION SYMMETRY 

Associated with the rotation symmetry, an additional 
Lagrangian is to be defined according to the gauge prin- 
ciple. It is important to observe from Sec. Ill A that, in 
the Lagrangian (19), the integrands of the last two in- 
tegrals are of the form 9t ( • ) • The action is defined by 
I = J / [At + 5t( • )] drd'^a. This property is regarded 
as the simplest representation of topology in the gauge 
theory ([16] ~ [18], [19]). In the context of rotational 
flows, it is known that the helicity (or Hopf invariant, 
[15]) describes non-trivial topology of vorticity fleld, i.e. 
mutual linking of vorticity lines. This is closely related 
with the Chern-Simons term (without third-order term) 
in the gauge theory. This term lives in one dimension 
lower than the original four space-time (a;^) of the action 
I because a topological term in the action is expressed 
in a form of total divergence (S^i^^) and characterizes 
topologically non-trivial structures of the gauge fleld. 

However, we learn here from the formulation of 
Sec. Ill A and look for a r-independent field directly. 



A. Lagrangian and helicity 

The r-independent field can be found immediately 
from Eq. (10). Taking the curl of this equation with 
respect to the coordinates (a, b, c), we obtain 

'^aXdrVa = dr{WaXVa) = 0, (30) 

where Va = {da,db,dc). Hence, one may write as Va x 

Va = ^a{a) [7]. 

The vector Va is a transformed form of the velocity 
V = {Xr,Yr, Zr) = {u,v,w) iuto the a-space. This is 
seen on the basis of a 1-form defined by 

yi = Vada + Vbdb + Vcdc (31) 
= udx + vdy + wdz. (32) 
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where Va = uxa + vya + wza^ Xa = dX/ da, u = Xt, etc. . 
Its differential dV^ gives a two- form il^ — dV^: 

= Q^db A dc + Ofcdc A da + l^cda A db 

= Uxdy A dz + ujydz A da; + Uzdx A dy, (33) 

where {Qa,^b,^c) = ^a, and V x t; = (wa:,^^, w^) = a; 
is the vorticity. Thus, it is seen that is the vorticity 
transformed to tlie a-spacc. The equation (30) is trans- 
formed into the r-derivative of the 2-form , Lq^ fl^ = 
(understood as the Lie derivative). 

Next, let us introduce a gauge-potential vector Aa = 
{Aa,Ab,Ac) in the o-space, and define its 1-form A^ by 
= Aada + Ab db + Acdc = A^ dx + Aydy + Az dz. 
Thus, it is proposed that a possible type of Lagrangian is 

Aa = - / {drA,,n,)d^a= [ {A,Ew[i^])d''x, 

where = dtuj + {v ■ V)a; — {uj ■ V)v + (V • v)u:. 

New results were deduced from this Lagrangian in [3] : 
(i) the velocity v includes a new rotational term, (ii) the 
vorticity equation is derived from the variation of A: 

Ew[u>] = dtui + {v ■ V)u> - {u> ■ V)v + (V ■ v)uj = 0, 

and (iii) we have non-vanishing helicity H, where 

rr f a3 f Ew[cm\A] 3 

H = u! • V d X = u) ■ d X, 

Jv Jv P 

B. Uniqueness of transformation 

Transformation from the Lagrangian a space to Eu- 
lerian x(a) space is determined locally by nine compo- 



nents of the matrix dx^ /daK However, in the previous 
solution considered in Sec.IIC, we had three relations 
(11) between v = {Xr,Yr, Zr) and {Va,Vb,Vc), and an- 
other three relations (13) between A — (Xt-,-, Ftt, -^tt) 
and (Aa, Ab, Ac)- Remaining three conditions are given 
by the equation (33) connecting lj = {uJx,^y,0Jz) and 
r2o(a) = {Q.a, f^b, r^c). For example, Q.a is determined by 

= ipbV dcz - dcy dbz) + u)y {dbZ dcX - dcZ dbx) 
+ ujz {dbX dcy - dcX dby). (34) 

There are three vectors (velocity, acceleration and vor- 
ticity) determined by evolution equations subject to ini- 
tial conditions in each space of x and a coordinates. 
Transformation relations of the three vectors suffice to 
determine the nine matrix elemets dx^ /da} locally. Thus, 
the transformation between the Lagrangian a space and 
Eulerian x(a) space is determined uniquely. [3] 



VI. SUMMARY AND DISCUSSION 

Following the scenario of the gauge principle of field 

theory, it is found that the variational principle of fluid 
motions can be reformulated successfully in terms of co- 
variant derivative and Lagrangians defined appropriately. 
The present variational formulation is self- consistent and 
comprehensively describes flows of an ideal fluid. 

In the improved formulation taking account of the ro- 
tational symmetry with additional equations of (33), the 
transformation relations of the three vectors (velocity, 
acceleration and vorticity) suffice to determine the nine 
matrix elements dx^ /da^ locally. Thus, the transforma- 
tion between the Lagrangian a space and Eulerian a; (a) 
space is determined uniquely. 
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